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We present a joint experimental and computational study of the effect of bacterial motion on micron-scale colloids contained in 
a two-dimensional suspension of B. subtilis. With respect to previous work using E. coli, here we introduce a novel experimental 
'set-up that allows us to realise a two-dimensional bacterial suspension insensitive to either evaporation or fluid flow. By analysing 
the mean square displacements of both bacteria and colloids, we confirm the existence of a crossover from super-diffusive 
behaviour at short time scales to normal diffusion at longer times. We also study the same two-dimensional system by means 
^^ of numerical simulations, using a suspension of self-propelled dumbbells or the Vicsek model, which has been previously used 
to study the dynamics of active particles. Our numerical results obtained with both models are in broad agreement with the 
experimental trends, but only the dumbbell simulations can match the experimental data quantitatively. The level of agreement 
we find suggest that steric interactions due to collisions are important players in determining collective motion of the bacterial 
,bath, and should complement hydrodynamic interactions in experiments. 
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1 Introduction 

Schools of fish, microbial colonies and bacterial suspensions 
provide intriguing examples of self-driven soft matter which is 
far from equilibrium even in steady state^i. Their dynamics 
display spectacular pattern formation and collective behaviour 
and are dominated by long-range correlations, as can be ob- 
served in flocks of birds or schools of fish—. Likewise, bac- 
.teria can swarm and concentrated films of swimming bacte- 
ria facilitate complicated self-driven flows which are several 
.times larger in magnitude than swimming velocities of indi- 
viduals^: this phenomenon is sometimes referred to as "bac- 
.terial turbulence". Understanding the mechanism underlying 
■ordering and patterning in suspensions of self-driven systems 
Jias now become a topical question in soft matter and biologi- 
cal physics. 

Besides providing an intriguing playground for ideas in 
non-equilibrium physics, self-propelled living matter may 
,well provide templates for the next generation of artificial 
■smart materials, which can be made, for instance, to swim 
^n a controlled way in aqueous solutions. 'First-principle' 
physics studies of bacterial suspensions and colonies have al- 
ready been successful at providing possible mechanisms for 
the patterning, for instance in chemotactic strains- . 

Dombrowski and coworkers have experimentally studied 
cell-driven motions in droplets of Bacillus subtilis—. Pedley 
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t Electronic Supplementary Information (ESI) available: Video of a 9% sus- 
pension of B. Subtilis seeded with a small percentage of fluorescently labelled 
cells. Epi-fluorescence and standard bright-field illumination are combined 
to visualise all the suspended cells, but only the fluorescent cells are used for 
the particle tracking. The field of view is about 225 ^um by 168^m in dimen- 
sion and the video was recorded at a rate of 10 frames per second. See DOI: 
10.1039/bOOOOOOx/ 



and Kessler observed the formation of convection patterns in 
sessile drops, due to Rayleigh-Taylor instabilities induced by 
the motion of cells towards the free surface, richer in oxygen 
concentration™. In pendant drops, the authors observed flow 
patterns which originated from hydrodynamics interactions. 

In a seminal work, Wu and Libchaber- measured the dif- 
fusion of colloidal tracers in a bacterial bath, and found that 
these behaved superdiffusively at early times, and exhibited a 
crossover to a diffusional regime at later times. The crossover 
between the two different regimes occurred at about 1 s, which 
is about 10 times smaller than the Brownian time over which 
a single colloidal sphere would diffuse a distance equal to its 
size. The proposed explanation for this phenomenology is 
that at early times colloids are pushed around by swimmers 
which leads to their superdiffusive behaviour Then, due to 
bacterial tumbling, many body hydrodynamics interactions or 
noise, both colloids' and bacteria's trajectories become essen- 
tially random, giving rise to diffusive dynamics. The effective 
diffusion coefficient associated with the latter regime is strik- 
ingly different in magnitude from the one which could be in- 
ferred from the Stokes-Einstein law for Brownian particles of 
the same size as bacteria. Still, the diffusion coefficient turned 
out to exhibit the same functional dependence on the colloid 
radius as in the Stokes-Einstein's expression^. 

In order to interpret their experimental data, Wu and Libch- 
aber proposed a Langevin-type model where the "collisional 
force" acting on the colloids due to the bacterial motion 
was exponentially correlated with a memory effect. In a 
subsequent work, Gregoire and coworkers- pointed out that 
Langevin equations would predict a non-realistic ballistic be- 
haviour at short time scales and suggested to use a model 
based on self-propelled particles, the Vicsek model, to study 
the collective dynamics of a swarm of bacteria^—. In their 



Vicsek model, each self-propelled particle moves in the di- 
rection of its orientation, which is aligned with the average 
orientation of all particles in its neighbourhood, with some 
added noise. Particles also interact via a repulsive poten- 
tial with passive beads suspended in the system. Using this 
model, Gregoire and coworkers qualitatively recovered the 
super-diffusive behaviour which was observed experimentally 
at short time scales. 

However, it has now been reahsed that the existence of the 
crossover behaviour found by Wu and Libchaber is a generic 
phenomenon, to be expected on very general grounds in mod- 
els of dynamics of a suspension of self-propelled particles, 
whether living or synthetic. For instance, more recently a 
number of researchers have proposed numerical models for a 
suspension of self-propelled particles at low Reynolds num- 
bers. Llopis and Pagonabarraga— performed lattice Boltz- 
mann simulations of self-propelled particles, finding that "hy- 
drodynamic noise", coming from the interaction between dif- 
ferent swimmers, is already enough to lead to diffusive be- 
haviour at large times, even in models not incorporating bac- 
terial tumbling, and showing a crossover from ballistic to dif- 
fusive motion. The latter result was confirmed by subsequent 
works, e.g. the Stokesian dynamics simulations presented in 
Ref.— , and the numerical work by Ortiz and coworkers—. The 
latter work proposed another minimal model for a microscopic 
swimmer, capable of capturing the leading order (dipolar) far- 
field hydrodynamics interactions between swimmers. Each 
swimmer is modelled as a rigid dumbbell made of two beads 
connected by a rigid body; the drag acts on both beads and 
the flagellum exerts a constant force on one of the beads in the 
direction of the dumbbells' unit vector and an equal and oppo- 
site force on the fluid. Within this model, a crossover between 
super-diffusive behaviour at short time-scales and diffusive at 
longer time-scales was found. At high concentration the au- 
thors observed large-scale coherent fluid motion. However in 
none of these works a suspension containing both passive and 
active particles was considered as in the original experiments 
in Ref.l. 

We note, in passing, that another interesting theoreti- 
cal framework within which to study suspensions of self- 
propelled particles is via a continuum theory, which has been 
successfully used to study microbial and bacterial fluids, for 
instance, in Ref.—"—. It would be interesting to explore what 
predictions these theories yield as regards the existence and 
behaviour of the crossover in the dynamics of colloidal tracers 
in active fluids - however here we do not address these ques- 
tions as these models work in the limit of high density whereas 
our experiments are with relatively dilute suspensions. 

Notwithstanding the significant amount of both experimen- 
tal and theoretical existing works, it still remains an open issue 
whether the early regime of the colloidal motion is ballistic or 
superdiffusive in nature. By considering a collection of non- 



interacting wild-type (run & tumble) bacteria, one might ar- 
gue that the early time regime should correspond to simple 
runs and be ballistic and the later times, after one or more 
tumbles, should be diffusive, as the bacteria are effectively 
performing a random walk at large length and time scales. 
However, most experiments show a superdiffusive behaviour 
at short time scales, with the mean square displacement be- 
having as a power law with an exponent a 1 .5 < a < 2^. The 
reasons behind this peculiar behaviour are still not understood. 
Furthermore, and importantly, we are not aware of a study 
which directly compares in detail experimental and simulation 
data. Such a comparison would obviously be useful, and for 
instance allow to more stringently and quantitatively test the 
dynamical predictions coming from models of self-propelled 
particles. E.g., are propulsion and collisions enough to quan- 
titatively match the data? Or is there a large effect of hydro- 
dynamic interactions coming from the (dipolar) forces which 
swimmers need to exert on the fluid in order to move? 

Our aim in this work is to attempt to answer these ques- 
tions by combining results from experiments and simulations. 
Experiments are performed on a quasi two-dimensional bacte- 
rial suspension of Bacillus subtilis and our numerical simula- 
tions adopt two independent ways to represent self-propelled 
particles. On the one hand, we simulate the bacterial bath as 
Brownian Dynamics of self-propelled dumbbells, interacting 
via hard repulsion with passive colloids. This allows us to ap- 
propriately consider the effect of steric interactions, as well as 
of the rod-like shape of bacteria. On the other hand, we use 
the two-dimensional vectorial Vicsek model^^i^ with the addi- 
tional assumption that both self-propelled and colloidal parti- 
cles have a finite size, represented by a short range hard-core 
repulsion, and interact via a soft repulsion with the passive 
colloids embedded in the suspension. 

Our emphasis is on the comparison between simulations 
and experimental data in rather dilute bacterial films: in this 
way we can assess the accuracy of the various models pro- 
posed in literature for self-propelled particles and bacterial 
suspensions, with respect to (i) the existence and location of 
the crossover between superdiffusive and diffusive behaviour 
and (ii) the time series of the mean-square displacements of 
both the colloidal and the active particles. In our simulations 
we disregard fluid-mediated interactions between two dumb- 
bells or self-propelled particles (although the orienting term 
coming from the Vicsek model may be thought as arising 
from hydrodynamic interactions among other things), and we 
will see that this is not crucial to either the existence of the 
crossover or its quantitative estimation. The mean-square dis- 
placements recorded from the simulations are only in semi- 
quantitative agreement with experiments, however the dis- 
crepancy may be due to other details which are not included 
in the simulations, such as the modelling of tumbling events 
in B. subtilis (see Section|2] and the discussion of our results). 



2 Experimental set-up 

Similar to the pioneering work of Wu and Libchaber— we ex- 
perimentally study a quasi two-dimensional bacterial suspen- 
sion. In our experiments we sandwich a dilute suspension of 
motile Bacillus subtilis in between two oxygen-plasma-treated 
glass cover-slides, creating thin samples (about 5—7 fjm in 
height) that can be imaged with high optical quality. Par- 
ticle tracking of individual particles is achieved by using a 
small fraction of fluorescently labelled tracers, either passive 
polystyrene beads or fluorescently labelled motile bacteria. 

B. subtilis is a peritrichously flagellated bacterium often 
found in soil, with a general morphology characterised by a 
rod-shaped body approximately 4 /jm in length and 1 /jm in 
diameter (just after cell division). Although its swimming 
motility is not well characterised, it is assumed to perform 
a "run-and-tumble" motion similar to that of the extensively 
studied E. coli—. periods of swimming in straight trajecto- 
ries ("runs") are separated by shorter "tumbling" events where 
cells randomly re-orientate their body. Runs typically last for 
~ 1 s during which B. subtilis swims with a speed of about 
30/jm/s. Although tumble events only last for few millisec- 
onds, they randomise the direction of the trajectories. As the 
combination of a large number of run-and-tumble events can 
be approximated by a random walk-like motion, it is natural 
to expect a diffusive-like motion over long time scales. 



staining techniques by non-specific binding—. Figure[T|shows 
the results obtained after staining the flagella. 
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Fig. 1 Fluorescently labelled 6. .sK&fifo (DS 19 19), where the dye 
(Alexa Fluor 488 C5) only stains the flagella. Top: Immobilised cell 
stuck to the cover-slip, highlighting the number of flagella and their 
shape. Bottom: a cell performing a run showing the formation of 
multiple flagella bundles. This peculiar flagellar arrangement during 
the run will be discussed in detail elsewhere^^l. 



2.1 Bacterial cultures and cell staining 

B. subtilis (DS1919) cells are prepared from frozen stocks, 
where an overnight culture is first prepared for 16 hours in 
Luria-Bertani (LB) broth at 30° Celsius and then displaced in 
an orbital incubator (with a shaker speed of 200 rpm). Next, 
B. subtilis cells are incubated with 100 /Jg/ml of spectinomycin 
(Sigma-Aldrich, S4014) in order to allow for fluorescent la- 
belling. We then harvest a 50/j1 sample from the overnight 
stock to inoculate 5 ml LB broth, which is incubated for 5-6 
hours in order to reach the mid-exponential growth phase of 
the bacteria. We then wash the sample once using a motility 
buffer containing 0.01 M of KPO4, 0.067 M of NaCl and 10""* 
M of EDTA (at pH 7.0) before re-suspending the bacteria in 
motility buffer at the desired dilution. Notice that after the 
washing step bacterial growth stops. 

To highlight individual swimming cells we adapt a flagella 
labelling technique developed by Blair et al.^. The B. subtilis 
cells are genetically modified to provide a high density of spe- 
cific binding sites for a fluorescent dye molecules along the 
length of its flagella filaments (for details of the genetic mod- 
ification and staining procedure see Ref.— ). The main advan- 
tage of this specific binding is that it allows to stain only the 
flagella filaments, thus avoiding light scattering from the cell 
body, unlike the one observed in E.coli, originated in standard 



The staining technique does not only allow to extract the 
average number of flagella on the B. subtilis cells used for 
our experiments, but also to estimate their length and char- 
acterise their general morphology. On average, our B. subtilis 
cells have (14.5 ± 3.5) flagella filaments per cell with an av- 
erage filament length of (8.5 ± 0.2) jum. Bright field imag- 
ing shows that the cell body is (6 ± 2) ^um long (consistent 
with a population of cells at different stages of growth) and an 
average run velocity is (30 ± 10) /jm/s (around 40/jm/s in a 
very diluted sample). It should be noted that labelling of the 
flagella provides much more information than required for the 
analysis presented in the following, for which much simpler 
labelling of the cell body would have been sufficient. Fluo- 
rescent imaging of flagella of swimming cells reveals novel 
insight about the swimming characteristics of both individual 
B. subtilis cells as well as their interaction in dense suspen- 
sions which will be presented elsewhere^. 

For the experiments on thin films described below, sam- 
ples consisted of a dilute suspension of (unlabelled) motile 
B. subtilis cells seeded with a small fraction (< 5%) of flu- 
orescent tracers: either fluorescently stained B. subtilis cells 
as tracers for active swimmers or yellow-green fluorescent 
polystyrene beads of 4/jm diameter (Invitrogen, F8859) act- 
ing as passive tracers. 



2.2 Sample preparation and imaging 

The experimental sample geometry is schematically repre- 
sented in Fig. |2] Two glass cover-slips are cleaned with ab- 
solute ethanol on an ultra fine microfiber cloth before being 
exposed to oxygen plasma and then dried with compressed 
air The main purpose of this treatment, which makes the cov- 
erslips hydrophilic, is to prevents bacterial cells and colloidal 
tracers from tethering to the glass. Next, we gently place 2 /j1 
droplet of the suspension at the centre of the large cover-slip 
and locate the smaller cover-slip on top of it, thus creating a 
thin layer of bacterial suspension. The two cover-slips are then 
sealed to prevent both evaporation and fluid flow (see figure|2]i. 
The sample thickness was measured by focussing the micro- 
scope on colloidal or bacterial particles occasionally stuck on 
the top and bottom surfaces or on the top and bottom of the 
drop edges. This led to a lower bound of ~ 5 /jm. In view 
of the size of the colloidal particles, or bacterial cell body we 
estimate the actual thickness of the bacterial suspension to be 
in the range of 5 — 7 /jm. 




Objective 



Fig. 2 Schematic diagram of the experimental configuration, where 
a suspension of bacteria is confined in a layer about 5iJm thick 
between two oxygen-plasma treated glass cover-slips (22 mm 
X 50 mm and 22 mm x 22 mm, Menzel-Glaser). 



supplementary material) of a suspension with about (|) = 0.09 
surface fraction (computed as the fraction between total bac- 
terial area and sample area, by assuming that a single bacteria 
covers an area of 4/jm^), where 5% of cells are fluorescently 
stained: in the image, it is even possible to identify the swim- 
ming modes of cells, where flagella dispersion indicates cell 
tumbling. From this image it is also apparant that there is an 
accumulation of bacteria near the drop edge. B. subtilis is an 
aerobic bacterium which moves towards the oxygen rich area 
near the egde, leading to an increase in bacterial density. How- 
ever, a steady-state is established quickly (~ 1 min) and by the 
time during data aquisition is started there is no more appar- 
ent change in bacteria density: the clustering affects a region 
of about 20 /jm close to the contact line which does not ham- 
per our measurements as our field of view extends ~ 100/jm 
away from the interface. 
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Fig. 3 Snapshot of a suspension of (f) = 0.09 of B. subtilis cells 
sandwiched between two oxygen-plasma-treated cover-slips, where 
5% of cells are fluorescently stained. Bacterial cell bodies appear 
black in the figure. 



The suspension near the edge of the suspension drop was 
imaged using an inverted microscope (TE2000-U, Nikon) with 
an oil immersion objective lens (Plan Flour 40 x OIL, NA 
1.3, Nikon). Images were captured using a cooled CCD 
camera (CoolSnapHQ2, Photometries) in conjunction with 
Metamorph software (Version: 7.5.5.0, Molecular Devices). 
For fluorescent imaging, a GFP filter cube (49002 ET-GFP, 
Chroma Technology Corp) was used and the excitation light 
was synchronised with the camera in order to minimise expo- 
sure. This fluorescence imaging provides high contrast images 
that allows for accurate particle tracking even if some back- 
ground bright-field illumination is still present (see figure |3]l. 
We recorded movies of the sample at 10 fps which were then 
analysed using a library of particle tracking routines written in 
IDI^. 

Fig.[3]shows a still frame from a movie (excerpt available in 



3 Computer simulations 

We perform computer simulations on a similar system, i.e. 
a suspension of self-propelled particles with surface frac- 
tion of 0.03 and 0.09, respectively, containing few colloidal 
particles (about 1% of the total number of particles simu- 
lated). In our work we choose two independent ways to rep- 
resent the self-propelled particles: a Brownian dynamics sim- 
ulation of self-propelled dumbbells, and a generalisation of 
the two-dimensional vectorial Vicsek model"-. Both simu- 
lations incorporate passive colloids, which interact with the 
self-propelled particles via purely repulsive interactions. Both 
self-propelled cells and colloids are chosen to have the same 
size ratio as in the experiments. 



In our experiments, colloids have a diameter of ae ~ 4 fxm. 
In a bulk sample they were found to diffuse at about D ~ 
0.1 jum^/s (in agreement with Stokes-Einstein law), whereas 
bacteria swim at a speed of vq ~ 30/jm/s. One may there- 
fore determine a dimensionless "Peclet" number as Pe = 
which with our numbers is then ^ 1 100. We note that the ef- 
fective diffusion for a 4 /jm colloid in a thin sample may be 
smaller than the bulk value, by a factor of ^ 0.5 — 0.6 consid- 
ering an effective thickness of 7 /urn—. Considering the error 
and the variations in the numbers above (as well as the varia- 
tion in bacterial cell size) it is more reasonable to consider a 
range of values for the Peclet number: comparison with ex- 
perimental data suggests that a range of Pe ^ 700 - 1400 is 
physically relevant. 



3.1 Self-propelled dumbbells 

The first avenue we propose to simulate a suspension of self- 
propelled particles with a small number of Brownian colloidal 
particles as tracers, is to use a system of "active" (in the sense 
of self-propelled) dumbbells. 

We equilibrate a quasi two-dimensional system of Nd=9QQ 
homomeric hard-dumbbells in an A^VT ensemble using cubic 
periodic boundary conditions in x-y (the slab is 1 particle di- 
ameter height in the z direction). Each dumbbell is made of 
two spheres with the same diameter a, and each sphere inter- 
acts with the other INj — 1 by means of a truncated and shifted 
Lennard- Jones potential 
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for r < rc {r = r,- — rj, with i,j = l,...,2Nd), whereas Vzj('") = 
for r > r^, where = 2'/^a and e defines the strength of 
the interaction. Having chosen such short range cutoff (that 
corresponds to the minimum of the interaction potential), al- 
lows us to eliminate the attractive contribution of the potential 
and make the spheres "effectively" hard. Two spheres within 
a dumbbell are also "glued" together by means of a stiff har- 
monic potential Vnir) ~ K(r — ro)', where ro = O is the equi- 
librium bond distance and /iT = ^ ^ = 1 Os/O" ( with the spring 
constant). 

Simulations are carried out using the open source 
LAMMPS Molecular Dynamics package (see on-line doc- 
umentation at http://lammps.sandia.gov/ and Ref.— ). The 
equations of motion governing the dumbbell motion in an im- 
plicit solvent (hence neglecting hydrodynamic interactions) 
are given by the (under-damped) Langevin equations acting 
on each sphere: 
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where m is the mass of a sphere, ^ the friction coefficient 
(^ = my with damping coefficient y), V the total conserva- 
tive potential acting on each particle {V = Vu + Vh) and Fr 
the force due to the solvent at temperature T - following the 
fluctuation-dissipation theorem which we expect to hold in the 
Fa = limit, we take this to be proportional to ^/kBT%R{t) 
with /?(r) a stationary Gaussian noise with zero mean and vari- 
ance {R{t)R{t')) — 5(r — t'). Self-propulsion is introduced via 
an extra force (Fa) acting on each sphere such that the force 
acting on the dumbbells has a modulus of 2|Fa| and the di- 
rection given by the dumbbell "director", which joins its front 
bead to its rear one (front and rear are chosen at the beginning 
of the simulation randomly). Such a director field corresponds 
then to the one used in polar nematic liquid crystals. 

We then add to the suspension of self-propelling dumbbells 
2%Nd of hard-sphere-like colloids, with <3c = 2a, and inter- 
acting with other constituents (passive or active) via the same 
truncated and shifted Lennard- Jones potential, where now a 
in eq[T]is half the sum of the hard core diameter of the two 
interacting particles. 

We start the simulations from an initial configuration where 
particles positions are random, while velocities come from a 
uniform random distribution at a given temperature of T = 
O.OOle/^fi. In our simulations, we sety=10T^' (lis the time 

unit, W with mass set to one), the time step 8f = 10^ T, 

and typically evolve the system for 10*' time steps. We choose 
the applied friction in order to simulate an effectively over- 
damped motion on the length scale of the particle size, as it 
is realistic for bacteria in our experiments. Parameters such 
as the temperature and the amplitude of the force Fq, are cho- 
sen in order to lead to a Peclet number comparable to the one 
found in experiments [Pe =1100 when the average bacteria 
velocity is 30/jm/s). More details on this model of active 
dumbbells and its properties will be presented elsewhere^. 

Our simulations of the self-propelled dumbbells are for 
some aspects similar to those performed by Ortiz and cowork- 
ers—, who did not have a specific data set as target for their 
work, and used a minimal model of swimmers able to capture 
the far-field hydrodynamics, focusing on collective dynamics 
in a confined suspensions. Our approach employs a purely 
Brownian dynamics simulations, therefore by comparing with 
the trends presented in figure 2 of Ref.— it is possible to single 
out the effect of hydrodynamic interactions. 

In what follows, we will indicate with a subscript b every 
variable that refers to the bacteria (or self-propelled particles), 
and with a subscript c every variable that refers to the col- 
loids. To quantitatively compare numerical results with ex- 
periments, we consider the diameter of the colloid Oc as unit 
length when measuring the mean square displacement of col- 
loids, o as unit length when measuring the mean square dis- 
placement of bacteria/self-propeUed particles, and the Brow- 



nian time to ^ (where D is the diffusion coefficient of the 
colloids) as the unit time for both colloids and bacteria. To 
make contact with experiments, the surface fraction (^i, is set 
to = 0.03 and (^i, ~ 0.09, respectively. 

Furthermore, we tune the value of the active particle veloc- 
ity vo such that the dimensionless number Pe ~ is kept 
constant and within the experimental range - we loosely refer 
to this as the "Peclet" number even though one should note 
that this ratio is defined with the diffusion of the colloidal 
particles and the velocity of the self-propelled ones. Such a 
choice is justified as D and vo give the natural diffusion and ve- 
locity scale in the model, as they are input to our algorithms, 
whereas for instance the diffusivity of bacteria can be found 
via numerical simulations. 

3.2 Vectorial Vicsek model 

As a simplified simulation to study the motion of bacte- 
ria as self-propelled particles, we also implement the two- 
dimensional Vicsek model. The Vicsek model was originally 
proposed to study flocking behaviours— as discussed in the 
original paper and later on e.g. in Ref.— . Here, we use the 
Vicsek model to study bacterial suspensions, in view of the 
facts that (i) bacteria can be considered as self-propelled par- 
ticles and (ii) they align with each other when they are close 
together, either via steric collisions or through hydrodynamic 
interactions^. The Vicsek model is an interesting example 
of an active non-equilibrium system, where activity leads (via 
self-propulsion) to the existence of an order-disorder transi- 
tion even in 2D, in violation of the Mermin-Wagner theorem 
which works for equilibrium systems. 

Here we use the so-called "vectorial" Vicsek model, whose 
dynamical rules are as follow s ^^i^" : the updated direction of 
the th particle is identified by an angle, 9, (f + Af), which 
is computed by (i) first adding noise to the direction of each 
particle within a given distance (ro) and (ii) averaging over 
these angles 



di{t + At) = arg 



,exp 



»^W+e«J,exp'^'(" 



(3) 



where arg denotes the argument of a complex number, and 
determines the direction of the velocity vector v,, is the 
number of neighbours within ro at time f including particle ; 
itself, ^,(f) is a delta correlated white noise (^/ ^ [ — 1, 1]) and 
e is the noise strength (e e [0, 1]). 

A relevant order parameter for the Vicsek model is the mean 
particle velocity ((p = \ Y.fjli '^i I- with vq the modulus of 
the speed and A^^ total number of Vicsek particles), as a func- 
tion of noise (e) and surface fraction {(^t>). 

Here we generalise the Vicsek model to incorporate a small 
number of colloidal particles (a similar simulation method was 



proposed in Ref.~, but no attempt was made there to get quan- 
titative predictions of mean square displacements). In our im- 
plementation, we introduce a finite size for both Vicsek parti- 
cle (which now have diameter Ob as opposed to being point- 
like) and passive colloids (with diameter 0^). 

In our algorithm, at every time step At each Vicsek particle 
evolves according to 

r,-.i(f+Af) = r,-.i(f)+v;,fc(f + Ar)Af (4) 
At At ^* 

j^i ,= 1 

where r, is the particle displacement, v,-.i is the particle ve- 
locity, = 3Kr[ah is the particle drag (r| is the viscosity) 
and Nc is the number of colloids; Uhc is a soft repulsive po- 
tential between the Vicsek particles and the colloids, whose 
gradient is constant and different from zero only when a par- 
ticle and a colloid overlap, whereas Ubb is a similar repul- 
sive potential between colliding Vicsek particles. We choose 
'^^+L,) = ^ 0.37. Furthermore, choosing voAf < ro 

guarantees that particles have enough time to interact with 
their neighbours. 

Passive colloids are modelled as two dimensional disks with 
diameter evolving via a similar over-damped Langevin dy- 
namics as in eqlH given by 

At ^* 

r,-,,(f + Af) = ri,,{t)- — Y,^Ubc{rij) (5) 
At 

- ^l^Vf/„.(r,7)+Fr 

V 1=1 

where r, ^ is the colloidal displacement, = 37tr|ac-, Fr is a 
random noise associated to the colloidal diffusion coefficient 
as in Eq|2] and Ucc is the repulsive potential between colliding 
colloids, (^^ ~ 0.37) and Ubc is the same as before. We 
simulate a two-dimensional square box with size L and peri- 
odic boundaries, setting Nh = 1000 and A^^ = 10 and choose 
Oc = 2a b- 

A crucial parameter within the Vicsek model is the amount 
of local ordering in the system, cp, which depends on a com- 
bination of the noise amplitude (e) the interaction range (ro) 
and Vo. When (p is large, Vicsek particles form large flocks 
and colloids move ballistically in the direction of the flock; at 
intermediate (p, Vicsek particles present low order and form 
short-lived flocks and colloids present a crossover (identified 
by a crossover time tc) from superdiffusive to diffusive be- 
haviour; for small (p, Vicsek particles are randomly oriented 
and colloids simply diffuse, albeit faster than purely Brown- 
ian particles. In our simulations, we choose an intermediate 
noise amplitude, £ = 0.55. This choice is motivated by the 
experimental observation that the suspension does not contain 



large flocks of bacteria, and that experiments show superdif- 
fusive behaviour at early times. Tuning ro also affects the be- 
haviour (at a given (p) and the position of the crossover time 
tc'. the smaller ro, the earlier tc occurs. In our simulations, 
we observe the best (qualitative, see below) match with the 
experimental trends with = 20^ (even though we also tried 
ro = 0( and '"o = 30^): our best results correspond to a value 
of the "Peclet" number (Pe = '■^) of the order of ^ 1400 to 
which we stick to in the following. 

4 Results 

To start our analysis, we show some typical snapshots of the 
bacterial suspension we study. Figure |4] shows some opti- 
cal microscope images showing both bacteria and fluorescent 
beads, at the two densities we have recorded our data for, 
(^b = 0.03 (top) and i^t> ~ 0.09 (bottom). It can be seen that 
in both cases there are significant inhomogeneities in the bac- 
terial density. Our experiments also suggest that there is col- 
lective dynamics of bacteria at both the densities studied (see 
Supplementary Movie 1 for a video of the dynamics of bac- 
teria at (^h — 0.09). Such a collective behaviour may be due 
to either collisions or hydrodynamic interactions between the 
bacteria. Our experiments and simulations however both sug- 
gest there is negligible nematic ordering in the sample at these 
densities. 

Figure|5]shows snapshots taken from our numerical simula- 
tions (where bacteria are either self-propelled dumbbells (red 
dumbbells in the bottom panel) or Vicsek particles (red arrows 
in the top panel). These snapshots also give evidence of a sig- 
nificantly inhomogeneous density profile which qualitatively 
matches the experimental one. 

In order to explore more systematically and quantitatively 
the dynamics of both colloidal particles and bacteria in exper- 
iments and simulations, we recorded time series of the mean 
square displacement (MSD) of both - the latter is defined via 
MSDh i;{t) = {[n{t) — n{Q)]^)ii,c, where the average runs over 
the total number of bacteria or colloids (the MSD computed 
from numerical simulations is then averaged over tens of real- 
isations). The mean square displacement depends on time as 
where a = 1 when particles move diffusively, 1 < a < 2 
when they move superdiffusively and a = 2 when they are in 
a ballistic regime. 

We first analyse the mobility and dynamics of the colloidal 
tracers. Figure |6] shows the mean square displacement of the 
colloids as a function of time at (^t> = 0.03 and = 0.09 
estimated from experiments and simulations where the self- 
propelling particles where either the self-propelled dumbbells 
(SP) or the vectorial Vicsek particles (in the inset). 

/,From Fig. |6] it is apparent that the MSD computed from 
the SP dumbbells matches the experiments better than the 




Fig. 4 (a) B. subtilis cells suspended between treated cover-slips at 
^l, = 0.03 (top) and ^i, = 0.09 (bottom), and 'X^m fluorescent beads. 



MSD from the Vicsek particles - the latter largely underes- 
timates the passive colloidal mobility. The quantitative dis- 
crepancy between our two simulations may be either due to 
the simplified steric representation of colloids and bacteria in 
the Vicsek model (which are both spherical) or to the presence 
of residual flocking, which may lead to unrealistic density in- 
homogeneities in the sample, possibly leading to an underes- 
timation of the collisions and of the interactions with the dif- 
fusing colloids. Nevertheless, we find that the trends of the 
curves obtained with the Vicsek simulations are qualitatively 
in agreement with the ones measured both experimentally and 
with the SP dumbbells simulations: when plotting the MSD of 
the colloids as a function of time, we observe a transition from 
a short time regime where colloids move in a superdiffusive 
way, to a long time regime where colloids move diffusively. 
Moreover, we observe that colloids move much faster (about 
2 order of magnitude) than their purely Brownian counterpart 
{MSD = 4tD shown as the green line in Fig.|6]l. 

The crossover times predicted by our simulations (both with 
dumbbells and Vicsek particles) together with those found ex- 




Fig. 5 Top: Bacteria particles simulated using the vectorial Vicsek 
model (red arrows) at surface fraction ^i, = 0.09. Passive beads are 
represented by green dots. Bottom: Bacteria simulated as 
self-propelled dumbbells (in red) and passive tracers (in green) at 
the same surface fraction. 



perimentally are summarised in Table [T] where we also in- 
clude the early and late time exponents characterising the 
MSD. 

Table 1 Results from experiments [Exp] and numerical simulations 
for colloids in self-propelling dumbbells [SP] (at different Pe 
number) and in Vicsek particles [Vic] at (^h = 0.03 and 0/, = 0.09: 
the crossover time tc is expressed in units of to = /Dc and a.s and 
«/ correspond to the exponent of for short and long time, 
respectively. The last column shows the value of the MSD scaled by 
CT at the crossover, MSD{tc)/al. Relative errors in the 
experimental and simulated estimates of the crossover times are 
about 10%. Typical errors on exponents do not exceed ±0.1. 
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Fig. 6 Beads MSD (in units of a^) as a function of time (scaled with 
fo = O^/Dc) at (^h = 0.03 (red curves) and ^i, = 0.09 (black curves). 
Open squares are experiments, filled diamonds results from 
simulations with SP-dumbbells. The thick green line shows the 
Stokes law prediction for the MSD for diluted Brownian colloidal 
particles in the absence of any swimmers (this agrees well with our 
data for isolated colloids in solutions, see above). Inset: Beads MSD 
from Vicsek simulations and experiments (both in units of Oc) as a 
function of time (scaled with fo) at the same ^i, (^i, = 0.03 (red 
curve) and (^i, = 0.09 (black curve)). Open squares are the 
experimental results (same as the ones in the figure). 



Note that in experiments the Peclet number is in the range 
[700 — 1400], so that simulations with SP dumbbells are done 
with Peclet number of 800 and 1 100 (simulations with the Vic- 
sek particles are shown for a Peclet number of 1400 which 
yields the best results). As shown in Table [T] there is a good 
quantitative agreement between the estimates of the crossover 
time, at both bacteria surface fractions, obtained from experi- 
ment and simulations. Even though the crossover time is well 
predicted, the simulations at all Peclet number slightly overes- 
timate the exponent associated with the superdiffusive regime. 
Quite possibly, the discrepancy may be due to our simplified 
description of the swimming mechanism of B. subtilis (see 
Discussion below). Our simulation data also suggest that the 
crossover time decreases when the Peclet number increases, 
at least in the range we considered. All the data presented in 
Table [T] both experimental and numerical, also show that the 
crossover time increases with the density of the system (al- 
though the experiments alone would not be conclusive as there 
is only a small change in the MSD exponent). In what fol- 
lows, we will run the simulations only for the system contain- 
ing dumbbells, and only at the Peclet number corresponding 
to the average bacteria velocity measured in the experiments 
(i.e.30^m/s). 

In order to investigate more in depth the dependence of the 
crossover time with density, we simulate several systems con- 
taining SP dumbbells and colloids, and systematically vary the 



dumbbell surface fraction, (^t>- Our results for the colloids are 
presented in Table|2](data of t^, a.s and a/ are obtained plotting 
the MSD of colloids). 



Table 2 Results for tc (in units of [a^/Dc]), ot.s and a/ for colloids 
immersed in suspensions with increasing (|)/,. Relative errors in the 
simulated estimates of the crossover times are about 10%. 
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The data presented in the table give evidence of an intrigu- 
ing non-monotonic relation between tc and (^t,. When the sys- 
tem is diluted (0.001 < < 0.03), colloids diffuse around 
and sometimes are pushed by the few swimming dumbbells in 
the system: in this range a higher density gives rise to more 
randomised collisions, hence the crossover to diffusion occurs 
earlier. However, for (^i, > 0.03 this trend reverses and an in- 
crease in density leads to an increase in the crossover time. 

To further understand the dynamics of our system, and to 
understand the mechanism leading to the non-monotonic trend 
of the crossover time with density, we set out to study more in 
detail the motion of the bacteria, or SP dumbbells. To this 
aim, we have measured, by means of both experiments and 
computer simulations, the mean square displacement of bac- 
teria/SP particles in bacterial bath at different surface fraction 
{<^b = 0.03 and (^b = 0.09) as well as in the dilute limit (in 
which bacteria are so dilute that they do not interact with each 
other). The results are shown in Fig.|7] 

As shown in the inset (for experiments only, simulations 
show the same trend), we also observe that bacteria (open 
squares) move much faster than the beads (thick continu- 
ous lines) at the same surface fraction ((|)i = 0.03 in red and 

= 0.09 in black). At short time scales, the data coming 
from simulations of the SP dumbbells agree well with the 
experimental ones. Qualitatively, the dumbbell simulations 
and the experiment also agree on the trend that the higher 
the density becomes, the slower the long-time diffusion of 
the bacteria/self -propelled particles is. This effect comes from 
crowding due to the steric effects of particles. Quantitatively, 
on the other hand, the SP particles move faster than the bacte- 
ria in the experiments. This may be due to additional hydrody- 
namic drag felt by the bacteria in the experiments, as they are 
swimming in a thin fluid layer (thickness ^5—7 /jm) sand- 
wiched between two cover-slips, so that wall effects may play 
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Fig. 7 MSD of bacteria/self-propelled particles (in units of a) as a 
function of time (scaled with the Brownian colloidal time to) at two 
bacterial surface fractions, d^/, = 0.03 (red curves) and ^i, = 0.09 
(black curves). The cyan curve represents the experimental results 
for a diluted system. Open squares are experimental results, filled 
diamonds are simulations with SP dumbbells. In the inset: beads 
MSD (from Fig[6]l and bacteria MSD (MSD is measured in fjm? and 
t in seconds). 



a role. 

The result that the MSD of bacteria at long time scales de- 
creases with increasing density, seems at first sight incompat- 
ible with the results reported in Ref.- which shows that, due 
to cooperative effects, dense suspensions of B. subtilis swim 
collectively /a^fer than single bacteria. One possible explana- 
tions is that our densities are too dilute in order for collective 
effects to play a role. But what is the mechanism leading to the 
onset of cooperativity? Does this arise from hydrodynamic in- 
teractions? Or is the steric exclusion enough for it to appear? 
In order to explore this possibility, we have simulated signif- 
icantly denser dumbbell suspensions which we were not able 
to study experimentally. The resulting MSD plots for the self- 
propelled particles are shown in Fig. [8] 

^From the MSD data in Fig. [8] one can estimate the (late 
time) diffusion coefficients of the dumbbells as a function of 
i^b- These are plotted in the inset of Fig. [8] where we ob- 
serve again a non-monotonic trend with this time of the 
SP particle diffusivity. First, as observed above, mobility and 
effective long time diffusion decrease with density, as it is in- 
creased from 0.03 to 0.06. After this, the long-time diffu- 
sion increases, and this enhancement becomes more and more 
dramatic with increasing density: for instance the MSD for 

= 0.15 is already an order of magnitude larger than the 
one at (^b = 0.03. Therefore our simulations suggest that even 
steric interactions alone can lead, through over-damped col- 
lision, to long range correlations and an enhanced mobility. 
What is particularly interesting is that this effect is already 
well developed at a relatively modest density of 0.15. This 




Fig. 8 MSD (in units of o) of SP dumbbells as a function of time 
(scaled with the Brownian colloidal time Iq) at various surface 
fractions. The inset shows the long-time diffusion (in units of /jirr/s) 
of the dumbbells as a function of (|)/, (the red squares are the 
experimental diffusion coefficients multiplied by a factor of 7 ~ this 
shows that in spite of the quantitative discrepancy the trend is the 
same). 



steric mediated cooperative motion is thus going to compete 
to the hydrodynamic cooperative motion, which is neglected 
in our Brownian dynamics simulations. 

Fig. [8] demonstrates that the interplay between steric in- 
teractions (crowding) and (purely collisional) collective be- 
haviour can give rise to a non-monotonic SP particle diffu- 
sivity with density. We believe that the same interplay is at the 
base of the non-monotonic trend of the crossover time for col- 
loidal particles as a function of To further consolidate this 
argument, we show in Table [3] the crossover time in the dy- 
namics of the SP particles: once more, this quantity behaves 
non-monotonically with (])/,. Table [3] shows the early and late 
MSD exponents, and a/, which we found from the MSD of 
SP dumbbell in suspensions at various (^h. 



Table 3 Results for tc (in units of [a^ /-Dd). OL.i and a/ for 
self-propelling dumbbells in suspensions with increasing ^i,. 
Relative errors in the simulated estimates of the crossover times are 
about 10%. 
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We observe that while in the most dilute case the dumbbells 
initially move ballistically, the early time behaviour crosses 
over to superdiffusive for increasing Note that for the most 
diluted system ((|)i = 0.001), the long-time exponent a/ indi- 
cates a superdiffusive behaviour, as the self-propelling dumb- 
bells so rarely collide with their neighbours - however for all 
higher surface fraction we reproduce the late time diffusive 
behaviour observed in experiments. 

5 Discussions and conclusions 

One of our main findings is that in both experiments and nu- 
merical simulations we observe that, at short times, passive 
colloids are dragged by bacteria and move in a super-diffusive, 
rather than ballistic, whereas at later times they slow down and 
diffuse, although their long-time diffusion coefficient is very 
different from, and larger than, their infinite dilution limit. The 
time at which the crossover occurs is a small fraction of the 
Brownian time which it takes for a particle to diffuse its own 
diameter, and its numerical value can be reproduced by our 
simulations with high accuracy. The presence of a cross-over 
from super-diffusive to diffusive behaviour is in agreement 
with Wu's experiments of a two-dimensional diluted bacte- 
rial suspension on a soap film—, where at short time scales the 
system turned out to be super-diffusive with (Ar^(f )) ^ r" and 
1 .5 < a < 2, slightly larger than in our experiments and closer 
to our numerical estimates^. 

Even though superficially our experimental work is similar 
to the one by Wu and Libchaber ^, there are several important 
differences. Firstly, the previous work used a different bac- 
terium, E. coli, which is smaller (typically 2x1 /jm) than the 
B. subtilis used here. This makes the bacteria in our experi- 
ment more similar in size to the passive tracer beads and, in 
combination with the thinner sample geometry, also ensures 
that the movement of our B. subtilis cells is almost strictly 
two-dimensional as cells cannot tumble in the direction per- 
pendicular to the cover-slips. Furthermore, as shown in Fig. 
1, the bundle arrangement of flagella while swimming is quite 
different from the one in E. coli, which will be discussed in 
a separate manuscript™. Secondly, the previous work em- 
ployed a soap film: soap is a chemical repellent which can dis- 
rupt swimming of bacteria and furthermore bacteria and beads 
can get stuck at the air-water interface of the film. Thin soap 
films are also prone to drift or rupture due to evaporation and 
capillary forces. In our experiments the motile bacteria sam- 
ple is instead suspended between two oxygen plasma treated 
cover-slips which makes them extremely robust against evap- 
oration and rupture. The oxygen plasma treatment also makes 
the glass surfaces very hydrophilic and thus prevents cells and 
beads tethering to the surface. The cover slips make it possible 
to use high resolution immersion objectives, improving parti- 



cle image analysis hence overall accuracy in the determination 
of trajectories and mean square displacements. Therefore, our 
setup is particularly well suited for quantitative comparison 
with computer simulations. However, the two nearby glass 
surfaces also introduce non-slip boundaries which may lead 
to hydrodynamic slow-down of the motion, as the plates are 
just about 5 — 7 /jm apart. This has not been taken into account 
in our simulations (or interpretation of the experiments). The 
experiments presented here rely on accurate tracking of the 
motion of individual particles. Even with the high contrast 
images obtained with fluorescence microscopy this is time- 
consuming and it is difficult to monitor enough tracers over 
a sufficiently long time interval to obtain good statistics. It 
has been recently demonstrated that better statistics can be 
achieved by means of novel experimental techniques such as 
Differential Dynamic Microscopy—. 

What is the superdiffusion-diffusion crossover physically 
due to? In their work, Wu and Libchaber concluded that its 
physical origin could be attributed to the collective dynam- 
ics of bacteria. But what ultimately leads to this collective 
behaviour? The latter might either be due to bacterial col- 
lisions, or to non-steric and hydrodynamic mediated interac- 
tions. While our dumbbell simulations faithfully consider the 
steric and geometrical properties of bacteria and colloids as 
well as their drag, long-range hydrodynamic interactions are 
neglected. Still, Fig. 8 clearly demonstrates that steric inter- 
actions alone (i.e. collisions) are able to lead to a collective 
behaviour, and also to the speed up of self-propelled particles 
at large densities, which was observed in concentrated bac- 
terial baths, and was commonly attributed to hydrodynamics 
alone. 

An important conclusion of our work is therefore that steric 
interactions via collisions play an important role, such that 
the collective behaviour is likely to come from a combina- 
tion of this effect and hydrodynamic interactions. An ac- 
curate model of hydrodynamic interactions in self-propelled 
particles is not a trivial issue, as accurate Stokesian dynam- 
ics simulations slow down dramatically with number of parti- 
cles^^^, and a more viable alternative might be the one sug- 
gested by Ref.— . Our algorithm employing Brownian dumb- 
bells is computationally less expensive and still leads to a 
good, semi-quantitative, agreement with experiments, which 
in many cases may be sufficient. 

Another interesting result from our analysis is that the Vic- 
sek model, although qualitatively leading to the existence of 
a superdiffusive-diffusive crossover and to an enhancement 
of the Brownian diffusion of colloidal tracers, shows a poor 
quantitative agreement with the experiments. To some ex- 
tent, this is quite surprising: while it is true that the Vicsek 
model does not resolve the rod-like shape of bacteria, it also 
does take into account the bacterial local alignment (which 
might be viewed as a rough way to mimic e.g. hydrodynamic 



effects) and has additional parameters which can control the 
local ordering of self-propelled particles - in dumbbells this 
only comes from geometry and steric effects. In view of the 
performance of the Vicsek model shown in Fig. 6, the compar- 
ison between dumbbell simulations and experiments becomes 
remarkable, given that there are essentially no free fitting pa- 
rameters in the simulations (except the Peclet number which 
experiments only constraint to be in a range of realistic val- 
ues). 

While the MSD of the colloidal particles coming from the 
simulations of the suspension of SP dumbbells and colloids 
can be quantitatively compared with the experimental ones, 
the MSD of the self-propelled dumbbells are only in order- 
of-magnitude agreement with the experimental bacterial ones. 
What is the reason for this behaviour? One possibility is 
that in our simulations we did not take into account bacterial 
tumbling, which may be thought to affect more the bacteria 
than the colloidal particles. However, it would not be triv- 
ial to modify the simulations to more faithfully reproduce the 
swimming mechanism of B. subtilis as this, to date, remains 
much less studied than that of E. coU - for instance, bacterial 
flagella bundle differently, as suggested by our experiments 
(Fig. 1 and Ref.~). Another possibility is that the neglect of 
hydrodynamic interactions is responsible for this quantitative 
discrepancy. Hydrodynamic interactions are thought to pre- 
dominantly come through the forces exerted by the bacteria 
on the fluid. These are not included in our treatment, whereas 
they may lead to an enhanced local flow. However a care- 
ful incorporation of the effect of such forces on the fluid may 
also not be trivial as it has recently been shown experimen- 
tally in-21 that the dipolar (stresslet) contribution is in reality 
difficult to measure, and masked by near-field quadrupolar in- 
teractions (and far-field monopole contributions due to the fact 
that swimmers in practice are not neutrally buoyant). 

Furthermore, for our experiments, due to the small thick- 
ness of the sample used experimentally, hydrodynamic inter- 
actions due to the flow field generated by the bacteria will be 
screened due to the presence of the nearby walls - this may 
partially account for the good level of quantitative agreement 
with our simulations. 

We should finally stress that it is quite possible that the 
effective density estimated in experiments is higher than the 
nominal one, due to inhomogeneities, which typically lead to 
a higher concentration at the edge of the droplet within which 
bacteria are confined. If the density were effectively larger, 
so would be the collective behaviour and the mean square dis- 
placement. 

In conclusion, we have performed a joint experimental and 
numerical study of the behaviour of colloidal tracers in a bac- 
terial fluid. With respect to previous simulation work, we have 
explicitly taken into consideration the finite size of both self- 
propelled particles (modelling bacteria), and of colloidal par- 



tides. The experimental work instead differs from previous 
work as we used B. subtilis films and a novel experimental 
setup. In agreement with earlier work, our experiments show 
that the effective diffusion coefficient of the colloidal parti- 
cles is about two orders of magnitude larger than the Stoke- 
sian value for 4/jm diameter spheres. By focusing on the 
traces of swimming bacteria and self-propelling particles, we 
show unambiguously that the early regime is superdiffusive, 
and we find that the crossover to diffusive regime occurs at 
a time which increases with bacterial cell density. We have 
compared experimental determinations of the crossover times 
with estimates from computer simulations, where the bacte- 
ria were modelled either via self-propelled (SP) dumbbells or 
via the (vectorial) Vicseck model. This is the first quantita- 
tive comparative study of this kind. The SP dumbbells sim- 
ulations reproduce the experimental trends qualitatively and 
the data semi-quantitatively even though they do not include 
long-range hydrodynamic interactions between the particles, 
thereby pointing to the importance of steric effects and colli- 
sions in setting up collective dynamics in self-propelled sus- 
pensions. The Vicsek model, while predicting qualitatively 
correct trends, leads to a rather poor quantitative agreement, 
perhaps due to the simplified steric description of the con- 
stituents of our system. Finally, our simulations show that the 
interplay of crowding and collective behaviour due to colli- 
sions can give rise to an intriguing non-monotonic behaviour 
of both the crossover time and the late time effective diffusiv- 
ity of self-propelled particles. It would be interesting to per- 
form similar joint theoretical and experimental analysis of sus- 
pensions of other synthetic, e.g. osmotically propelled^—, 
microscopic swimmers. 
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